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ABSTRACT 

The  author’s  conjecture  concerning  the  knot  sequence  whose  associated  B-spline  se¬ 
quence  has  maximum  max-norm  condition  number  is  disproved.  Related  condition  num¬ 
bers  axe  explored  and  the  corresponding  conjecture  concerning  the  ‘worst’  knot  sequence 
for  them  is  further  supported  by  numerical  results. 
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The  exact  condition  of  the  B-spline  basis  may  be  hard  to  determine 

Carl  de  Boor 

At  the  end  of  a  long  discussion  of  the  linear  functionals  which  vanish  at  all  B-splines 
but  one  in  [B2],  I  conjectured  that  D *i00,  the  worst  possible  condition  with  respect  to  the 
max-norm  of  a  B-spline  basis  of  order  k ,  occurs  when  the  knots  have  high  multiplicity.  I 
went  further  than  that  on  p.155  of  [B3],  where  I  displayed  supposed  values  of  Dfc.oo  based 
on  this  conjecture.  The  conjecture  was  based  in  good  part  on  detailed  calculations  of  a 
closely  related  problem  in  [B x ] ,  on  a  calculation  of  the  number  Dk  which  provides  a  bound 
for  the  worst  B-spline  condition  with  respect  to  any  p-norm,  and  on  some  calculations 
of  the  max-norm  condition  itself.  In  particular,  I  wrote  :  “As  with  the  earlier  reported 
calculations  of  Dk ,  it  appears  from  these  calculations  that”  the  worst  condition  “is  taken 
on  at  the  ‘middle’  vertex  of  the  simplex”  of  knot  sequences  over  which  the  maximization 
takes  place.  “This  would  mean  that 

=  IKAWeOr'lU 

with  r  :=  (r,)*fc  given  by 

0  =  TX  =  .  .  .  =  Tfc,  Tfc+1  =  .  .  .  =  r2jfc  =  1 

and  0  =  Qx  <...<£*=  1  the  extrema  of  the  Chebyshev  polynomial  of  degree  k  —  1 
for  [0, 1],  This  gives  the  following  values  for  Dk,oo  '■  . .  In  other  words,  I  conjectured 

that  the  worst  max-norm  condition  occurs  for  a  knot  sequence  without  interior  knots  and, 
assuming  this  to  be  true,  computed  and  displayed  this  condition  number  for  the  first  few 
values  of  the  order  k  as  the  value  of  Dk,0 o- 

Note  that  ||  (Nj.k,r(Qi))  ||oo  =  IMIoo  with  CT  :=  N i,k,rd{j )  the  unique  spline  sat¬ 

isfying  Cr{Qj )  =  {  —  j  —  1,  •  • . ,  k.  This  implies  that  Cr  is  the  Chebyshev  polynomial 
of  degree  k  —  1  for  the  interval  [0, 1]  and  the  numbers  computed  and  displayed  as  Dk,oo 
in  [B2]  and  [B3]  are  therefore  the  absolutely  largest  coefficients  in  the  expansion  of  the 
Chebyshev  polynomial  as  a  linear  combination  of  the  B-splines  N jtk,r-  Because  of  the  spe¬ 
cial  nature  of  the  knot  sequence  r,  these  B-splines  reduce,  on  the  interval  [0,1]  of  interest, 
to  the  polynomials  in  the  Bernstein  form.  This  led  Lyche  [L]  to  the  observation  that  there 
was  no  need  for  numerical  calculations  since  the  Bernstein  form  for  the  Chebyshev  poly¬ 
nomial  could  be  written  out  explicitly  and  a  simple  expression  for  its  absolutely  largest 
coefficient  could  be  provided.  Because  of  this  connection,  I  shall  refer  to  the  knot  sequence 
without  interior  knots  more  briefly  as  the  Bernstein  knots.  The  explicit  formula  allowed 
Lyche  [L]  to  verify  my  conjecture  that  this  condition  number  grows  like  2k . 

Since  then,  there  have  been  several  attempts  at  verifying  the  conjecture  that  the  worst 
max-norm  condition  is  had  by  the  Bernstein  knots.  It  is  therefore  important  to  point  out 
with  the  aid  of  specific  examples  that  the  conjecture  is  incorrect  in  general.  In  contrast, 
more  detailed  calculations  concerning  the  related  number 

Dk  ■=  sup  sup  sup  \a(j)\\Ij\/  /  | 

t  j  a  Jij 


1 


(with  Ij  [tj,tj+k]  the  support  of  Nj)  have  so  far  failed  to  shake  the  corresponding 
conjecture  that  Dk  is  attained  by  the  Bernstein  knots. 


Condition  number  defined 


It  is  convenient  to  define  the  condition  number  cond  of  the  basis  (<fi)  of  a  normed 
linear  space  S  as  the  number 


cond  :=  sup 


E.-y»'a(OII 


a  . 


a 


S^P  ij  tt-'  /  ■> 


For,  assuming  the  basis  (yq)  so  normalized  that  the  first  supremum  is  1,  this  gives  the 
equality 

cond  =  sup  ||  Ay  || , 
i 

with 

A  j  :  y ^<Pia(i)  a(j) 

i 

the  jth  coordinate  functional  for  the  basis 

Let  t  :=  ( t{ )  be  a  knot  sequence  for  splines  of  order  fc,  i.e.,  ti  <  all  i,  and  let  ( Nt ) 
be  the  corresponding  (normalized)  B-spline  basis  for  the  spline  space  S  :=  Sk,t  (see,  e.g., 
[B3],  for  relevant  definitions  and  details).  The  IV)  are  nonnegative  and  sum  to  1,  hence 


sup 

a 


II  E.^°(0IU 

llal|oo 


1, 


where  here  and  below  we  take 


in  case  t  is  finite.  Denote  by  A* 

condj^t 

its  condition  number. 


Il/lloo  :=  sup  |/(t)| 

tk<t<tn+1 

=  Ai,t  the  zth  coordinate  functional  for  this  basis  and  by 


=  sup 


11° 


a  II 


=  sup  ||Ai,t| 


A  counterexample 

Let  f  be  the  even  piecewise  cubic  on  [-1,1]  with  just  one  breakpoint,  at  0,  given  by 


jT3(l+(l-or)(*-l)),  x>0; 
nx)-\  x  <0, 
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with  T3  =  4()3  —  3()  the  cubic  Chebyshev  polynomial  and  a  :=  —1/2  its  negative  extreme 
point  (see  Figure  1). 


Figure  1.  A  cubic  Chebyshev  spline  with  one  knot  constructed  from  the 
cubic  Chebyshev  polynomial 

Since  .0/(0+)  =  0,  /  is  in  C2,  i.e.,  a  cubic  spline  with  a  simple  knot,  at  0.  One  readily 
computes  its  cubic  B-spline  coefficients  (for  the  knot  sequence  t:=(-l,  -1,  -1,  -1,  0,  1,  1,  1, 
1))  to  be  (1,  -7/2,  11/2,  -7/2,  1).  Since  ||/|joo  =  1,  this  implies  that  cond4>t  >  5.5,  while 
the  cubic  Bernstein-knots  condition  number  is  5  (see  [B2],  [L]). 


The  number 


A  lemma 


Dk,  00  :=  sup  sup  l/dist00)(ti.+1)t,.+t_l](iVi,  span  (IVj)^.) 

was  introduced  in  [B2]  as  a  convenient  upper  bound  for  the  worst  B-spline  condition 
number 

condfc  :=  supcond^^. 

t 

The  following  Lemma  shows  that  the  two  numbers  are  equal,  hence  that  cond*  =  Ojt;00 
can.  be  determined  by  local  means. 

Lemma  1. 

sup  sup  ||Ai||  =  sup  !|Afc_i,sj|/, 

t  i  s 

where  s  is  any  knot  sequence  of  the  particular  type 

Sj  =  .  .  .  =  5*  =  —  1  <  Sjfc+l  <  .  .  .  <  S2k- 3  <  1  =  •S2fc-2  =  .  .  .  =  S3fc_3 

and  || A|| 7  :=  sup f  |A/|/||/|]/  with  ||/|j/  the  max-norm  on  I  :=  [-1,  lj. 

Proof.  It  is  sufficient  to  prove  that,  for  any  t  and  any  i , 

||Aj||  <  D  :=  sup  1 1 A — ! ,* 1 1 / . 
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Since  A,()°  =  1  for  any  i,  we  have  min,  |j At ||  >  1  and,  in  particular,  D  >  1.  If  f,+ j  =  ti+ Jt_1) 
and  without  loss  of  generality,  <  U+i,  then  A,/  =  f(t{+ 1),  hence  ||Aj|]  =  1  <  D,  and  we 
are  done  in  this  case. 

In  the  contrary  case,  i*+i  <  £,•+*_ 1,  hence,  after  a  linear  change  of  the  independent 
variable,  we  may  assume  that  =  — l,ti+fc_i  =  1.  Now  let  t'  be  the  knot  sequence 
obtained  from  t  by  inserting  both  -1  and  1  enough  times  to  increase  their  multiplicity  to 
k  —  1  and  let  i'  be  such  that  =  ti+j  for  j  =  1,  ...,k  —  1.  Then,  with  A'  the  jth 

coordinate  functional  for  the  basis  ( Niyk,t ')  of  the  refined  spline  space  of  the  same  order, 

l|A,||<||A',||, 

since  the  (now  standard)  formula  (cf.,  e.g.,  p.116  of  [B3]) 

r<fc 


with  £,•  €]<,-,  <,-+jfc[  and 


</>  :=  (ti+ 1  -  -  • )/(k  -  1)!, 

shows  that  A,/  only  involves  the  knots  t,+y,  j  =  1, . . . ,  k  —  1,  hence  A ;/  =  A',  /  for  all  /  G  S. 
This  finishes  the  proof  since 

|| A',  ||  =  sup|A',/(/||/[(<x>  <  sup  |A(//|/||/j|/  <  D. 

!  f 


4k 


Computation  of  D jtj00 

According  to  Lemma  1,  D) t)00  is  the  maximum  of  the  function 

d  :  [—1,  l]fc  3  — +  IR  :  a  1— ►  ||Ajfc_ijS||/ 

with  (sjfe+.OfTi3  the  sequence  of  ‘interior’  knots  of  the  knot  sequence  s  =  (s,)3fe_3  in  [-1,1] 
obtained  from  a  by  ordering.  The  failed  conjecture  amounts  to  the  statement  that  the 
maximum  is  taken  on  at  the  ‘middle’  vertex  of  the  domain  of  d. 

For  k  =  3,  there  are  no  interior  knots  and,  correspondingly,  Dit00  =  3,  the  condition 
number  of  the  Bernstein-knot  B-splines. 

For  k  =  4,  there  is  just  one  interior  knot,  hence  the  calculation  of  Dkl0o  amounts 
to  the  maximization  of  the  function  d(£)  as  £  traverses  the  interval  [-1,1].  A  drawing  of 
this  function  is  available  in  Figure  2;  it  is  the  hindmost  curve.  This  shows  that  d  has  a 
local  minimum  at  £  =  0  (necessarily  a  critical  point  because  of  symmetry),  and  that  the 
maximum  (and,  at  least  numerically,  a  good  estimate  for  Dii00)  is  5.5680...  which  occurs 
when  £  ~  ±.472. 
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Figure  2.  Condition  number  of  cubic  B-splines  with  two  interior  knots. 

Sections  axe  shown  corresponding  to  one  knot  fixed  while  the 
other  traverses  [-1,1].  As  the  ‘fixed’  knot  traverses  [-1,1],  the 
corresponding  section  is  increasingly  offset  to  provide  ‘insight’ 
into  the  surface. 


It  is,  in  some  sense,  not  too  surprising  that,  for  k  =  4,  the  maximum  occurs  in 
the  interior  rather  than  at  a  vertex,  since,  after  all,  £  =  0  is  necessarily  a  critical  point, 
by  symmetry,  and  there  are,  correspondingly,  two  ‘middle’  vertices.  It  is  much  more 
discouraging  that,  for  k  =  5,  the  maximum  is  also  taken  on  at  an  interior  point,  for,  in 
this  case,  there  is  only  one  ‘middle’  vertex.  Figure  3  shows  d  as  a  function  of  the  two 
interior  knots.  According  to  [B2],  [L],  the  max-norm  condition  in  the  quartic  case  for  the 
Bernstein  knots  is  11  2/3.  But  one  computes  in  this  case  that  D 5<<x  ~  12.088  and  this 
occurs  when  the  two  interior  knots,  both  simple,  are  at  the  symmetric  points  ~  ±.89. 

The  sharp  drop  toward  the  boundary  values  is  an  indication  of  the  general  situation. 
Numerical  experimentation  for  k  <  8  seems  to  indicate  that,  for  k  >  3,  the  maximum 
occurs  at  an  s  close  to  but  not  at  a  vertex,  with  d  raising  sharply  initially  as  one  moves 
away  from  the  boundary.  For  odd  k,  the  maximum  seems  to  occur  near  the  ‘middle’vertex. 
For  even  k ,  it  occurs  at  a  point  (two  points  for  small  k )  near  what  passes  for  the  ‘middle1 
vertex  in  that  case,  i.e. ,  with  both  0  and  1  knots  of  the  same  multiplicity  and  1/2  a  simple 
knot. 

For  the  evaluation  of  HAjHoo,  consider  the  ‘Chebyshev  spline’  Ca  for  the  knot  sequence 
s,  i.e.,  C»  €  S  :=  Sk,a,  of  max-norm  1,  and  maximally  alternating ,  i.e.,  there  is  an  increasing 
sequence  ( Qj )"  (with  n  :=  dim  S)  so  that  C3(gj )  =  (— )""■>,  all  j.  (It  can  be  shown  that 
such  Ca  exists,  and  uniquely  so  (see,  e.g.,  [M].)  Let  c  be  the  sequence  of  its  B-spline 
coefficients.  This  sequence  necessarily  strictly  alternates  in  sign  at  least  n  —  1  times,  hence 
all  c(j)  are  nonzero. 
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Figure  3. 


Condition  number  of  quaxtic  B-splines  with  two  interior  knots 


This  implies  that,  for  each  j,  Cj  :=  Ca/c(j )  is  well  defined  and  in  NjtS  +  span (iVj)^, 
therefore  necessarily  the  error  in  the  best  uniform  approximation  to  N}  from  span (Ni)&j  = 

kerAJ  S,  hence  an  extremal  for  A JiS,  and  therefore  HA^doo  =  mcA  OO  -  This 

reduces  the  evaluation  of  the  function  d  to  be  maximized  to  the  numerical  construction  of 
the  Chebyshev  spline,  as  is  done  in  the  following  MATLAB  (cf.  [MBLK])  script. 

function  [sp , rho , a , iter] =chebmk (t ,  k , rho) 

X 

X  [sp ,  rho ,  a ,  iter]  =chebmk(t ,  k  [ ,  rho] ) 

X 

X  returns  the  Chebyahev  spline  for  the  given  knot  sequence  t_l,  ....  tjn+k, 

X  as  well  as  the  sequence  rho  of  its  alternating  points  and  the  sequence 
X  a  of  its  B-spline  coefficients.  On  input,  rho  is  assumed  to  contain  a 
X  reasonable  first  guess.  If  missing,  the  knot  averages  are  used. 

X 

X  By  definition,  the  Chebyshev  spline  is  the  unique  linear  combination 
X  of  the  B-splines  for  the  knot  sequence  t  which  has  norm  1  on  [t_k,t_n+l] 

X  and  takes  on  the  values  1  and  -1  alternatingly  the  maximum  possible  number, 

X  i.e.,  n  times,  and  is  positive  near  t_n+l  . 

npk=length(t) ;n=npk-k; 

t=[t(k)*ones(l,k)  t(k+l:n)  t(n+l)*ones(l,k)] ; 

X  here  I  omitted  statements  which  would  initialize  rho  as  the  average  of 
X  k-1  neighboring  knots  in  case  the  initial  guess  provided  is  inadequate. 

rho(l)=t(k) ;rho(n)*t(n+l) ;X  the  first  and  last  rho  are  the  endpoints  of  the 

X  interval. 

trho=rho(2:n-l) ;  X  only  the  interior  rho  will  be  iterated  on. 

y^onesCrho);  X  set  up  the  oscillating  data  to  be  matched  by  . . . 

y(n-l: -2: l)=-y(n-l :-2: 1) ;  X  ...  the  Chebyshev  spline. 
change=l;tsize=rho(n)-rho(l) ;  X  set  up  convergence  control. 

iter=0; 
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while  (change>l.e-8)*(iter<8) ; 

sp=spint(t,rho,y) ;  %  compute  the  epline  with  knot  sequence  t  which  takes 

X  on  the  value  y(j)  at  rho(j)  ,  j=l . n. 

dsp=spder(sp) ;  X  construct  the  first  derivative  of  this  spline,.. 

drho=spval(dsp,trho) ;  X  ...and  evaluate  it  at  the  interior  rho  . 
ddrho=spval(spder (dsp) ,trho) ;  X  also  evaluate  the  second  derivative  of  that 

%  spline  at  the  interior  rho  . 
drho=-drho./ddrho;  X  compute  the  Newton  step  ... 

trho=trho+drho;  X  ...  and  add  it  to  the  current  interior  rho  . 

X  prevent  modified  rho  from  violating  the  expected  interlacing  by  pulling 

X  back  on  the  proposed  Newton  step  if  necessary: 

count=0; 

while  Cany(trho<t(3:n)) |any(trho>t(k+l:n+k-2)) JanyCdiff (trho)<=0)) , 
drho=drho/2 ; trho=trho-drho ; 

count=count+i ; if  (count>20) , error( ’no  convergence') , end 

end 

change-max(abs(drho))/tsize  X  compute  relative  size  of  the  step  taken. 
rho(2:n-l)=trho;  X  update  rho  . 

iter=iter+l ; 

end 

[duffimy,a]=spunmk(sp) ;  X  recover  the  B-spline  coefficients  a  of  the  Chebyshev 

X  spline. 


The  calculations  become  quite  delicate  with  increasing  k  and  increasing  nonuniformity 
of  the  knot  sequence.  I  have  not  found  a  certain  rule  for  choosing  a  satisfactory  first  guess, 
but  have  very  often  succeeded  with  the  aid  of  continuation.  For  example,  if  the  Chebyshev 
spline  for  the  same  (interior)  knots  but  of  one  order  lower  is  already  available,  then  the 
midpoints  between  its  neighboring  extreme  points  often  provide  good  first  guesses  for  the 
interior  extreme  points  of  the  Chebyshev  spline  to  be  computed. 

Note  that  the  ‘Chebyshev  spline’  used  here  is  in  general  different  from  the  ‘Chebyshev- 
Euler  spline’  used  in  Schoenberg  and  Cavaretta’s  solution  [SC]  of  the  Landau  problem  on 
the  halfline,  and  which  also  appears  prominently  in  Tikhomirov’s  work  (cf.  [T]).  The  latter 
might  be  called  ‘perfect  Chebyshev  splines’  since  they  are  Chebyshev  splines  whose  highest 
nontrivial  derivative  is  absolutely  constant,  a  feat  achieved  only  by  an  appropriate  choice 
of  knots.  The  more  general  Chebyshev  splines  of  interest  here  have  most  recently  appeared 
in  Demko’s  [D]  nice  proof  of  the  existence  of  ‘good’  interpolation  points  for  arbitrary  knot 
sequences  and,  almost  simultaneously,  in  Mprken  [M],  a  reference  of  which  I  became  aware 
only  recently.  Mprken  devotes  an  entire  chapter  to  the  Chebyshev  spline  (which  he  calls, 
perhaps  more  helpfully,  the  fequioscillating  spline’),  proving  its  uniqueness  by  a  detailed 
study  of  the  sign  jtructure.  I  note  that  uniqueness  can  also  be  deduced  from  the  fact 
(mentioned  earlier)  that  Cy  is  an  extremal  for  Ay. 


The  1-norm  condition  number 


When  the  norm  on  S  =  5jt,t  is  the  1-norm, 


it  is  preferable  to  use  also  the  1-norm  instead  of  the  max-norm  for  the  B-spline  coefficients 
and  to  use  a  different  normalization  for  the  B-splines,  too.  Precisely,  define 

cond£it:=  ||*||||*-1||, 
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with 

$  :  — >  Sk  i  C  Li[<jk,in+1]  :  a  Mja(j) 

j 


and 

Mi  :=  := 

Since  ||A/j||i  <  1  for  all  j,  we  have 


tj+k  tj 


a  i 


hence 

:=  supcondj.  t  =  sup  ||4>_1||. 

t  ’  t 

It  is  worthwhile  to  point  out  that,  by  duality,  this  number  coincides  with  the  Favard 
constant  [Bx] 

,  inf{||D‘/IU  :  /  €  L <£>,/  =  /„  on  t} 

A (k )  :=  sup - — - 7 — - 

f0,t  maXjfc!|[<j,...,ti+jk]/o| 

which  measures  how  small  one  can  make  the  fcth  derivative  of  an  interpolating  function 
(relative  to  the  kth.  divided  differences  [t,, . . .  ,  t,+fc]/o  of  the  given  data).  This  fact  has  also 
been  found  by  Otto  [0],  by  rather  different  means. 

Lemma  2.  Dk,  1  =  K(k). 

Proof.  Since  k\[tt, . . .  ,<,•+*]/  =  f  M{Dkf,  we  can  write  K(k)  also  as 


K(k)  =  sup  sup 

t  So€Loo 


infdbiU  :  f  M}g  =  f  Mjg0,  all  j} 
maXj  \fMjgol 


This  shows  that 

K{k)  =  sup  j|F||, 

t 

with 

F  :  £oo  — *  S*  :  a  1-*  ^ 

and  nx  :=  fii,k,t  '•=  — 1  ^i,k,t  the  linear  functionals  dual  to  the  Mj's.  But  this  implies 

that  F  =  ($*)-1  and,  in  particular,  ||Fj)  =  jj$-1|j.  ♦ 

It  follows  that  the  calculations  of  A' (it)  in  [Bj]  are  pertinent  for  the  calculation  of 
Dk,i  ■  These  calculations  are  based  on  the  fact  that 


K(k)  <  1  +  2(fc  —  1) sup  ||A||, 

S 
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Figure  4.  ||A||  for  k  =  5  as  a  function  of  two  interior  knots 


with  0  =  s i  =  . . .  =  Sk  <  Sk+i  <  . . .  <  1  =  S2*-i  =  •  •  •  =  $3k—2  and  A  the  linear  functional 
on  Sk, a  C  Li[0, 1]  which  carries  Yhj  Mja(j)  to  see  [®i]  f°r  details. 

It  is  possible  to  compute  ||A||  as  a  function  of  (s,t+y)i_2  by  constructing  the  unique 
absolutely  constant  step  function  h  on  [0,1]  with  2k  —  2  jumps  for  which  A /  =  f  hf  for  all 
/  £  S.  The  calculations  are  almost  identical  to  those  reported  in  the  final  section.  They 
show  that,  for  small  k,  the  supremum  is  achieved  at  one  of  the  vertices  of  the  domain  over 
which  the  supremum  is  taken,  i.e.,  when  there  are  no  interior  knots.  This  is  illustrated  in 
Figure  4  for  k  =  5. 


The  p-norm  condition  number 

Finally,  consider  the  condition  number  of  the  B-spline  basis  when  the  norm  on  5  is 
the  p-norm, 

ll/ll,  :=  (  /'"+‘  \m,dl)'h. 

Jtk 


9 


It  is  shown  in  [B0]  (see  also  [B2])  that 


n-k'\\El,ra\\P  <  il  X>,aO')llp  <  HE'^a ||, , 

) 

Die  :=  sup  sup  A:  |I  jUi  j| (t) ,  (£aa)(j)  :=  (|/7|/fc)“a(j),  all  j, 

t  i 

and 

11  n(0  l^/l 
IlMir  ;  :=  SUPT-T^- 

/  J/,  I/I 

In  particular,  Dk  is  an  upper  bound  for  both  Dk,oo  and  Dk,  1,  while  also  kDk,  1  >  £>/;. 
Since  t  —  ^-±^— ^-A j,t,  it  follows  that 

Dk  =  sup||Afci.||, 

S 

with  =  .  .  .  =  Sfe  =  0  <  S/fc-,-1  <  .  .  .  <  s2/t-l  <  1  =  $2*  =  •  •  •  =  53*-l- 
Explicitly, 

in  lAot/l 

=  sup  sup  I (  sup  77 ■ 

t  i  f  J/{  I J  I 

Hence,  after  a  linear  change  of  variables  which  carries  the  typical  knot  interval  I,  =  tl+k] 
to  the  unit  interval, 

Dk  =  sup  sup —j—, 

8  /  Jo  I/I 

with  0  =  Sfc  <  s,t+i  <  . . .  <  S2k  =  1-  Since  this  involves  the  norm  of 

A  :=  AfciS 

on  S  :=  Sk, a  n  Li[0, 1],  the  knots  Sj  for  j  <  k  or  j  >  2k  are  immaterial;  I  take  them  to  be 
0  or  1  respectively.  The  remaining  knots  lie  in  [0,1]. 

Let  n  dim  5.  According  to  [B2],  ]J  A  J]  is  computable  as  the  absolute  height  ||^||oo  of 
the  unique  absolutely  constant  step  function  h  on  [0,1]  with  n  steps  which  represents  A  in 
the  sense  that 

A  /  =  /  hf  for  all  /  €  S. 

If,  more  precisely,  0  =  Sk+r  <  Sfc+r+i,  then  ||^j|oo  equals  the  norm  of  A  =  A^t  on  5  = 
Sk, t  H  Li [0, 1],  with 

t\  —  ...  —  t k  =  0  <  tk+ 1  ^  ^  in  1  =  tn+ 1  =  •  •  •  =  tn  +  * 

and 

£=  k  —  r. 
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The  following  MATLAB  script  returns  this  step  function  h  for  given  t  and  given  intt  := 

{tk+r+  1  i  •  •  •  >  in)- 

function  [beta , tau , iter] =stepmk(lef t , intt ,k , tau) 

X 

X  [beta, tau , iter] =stepmk (left , intt ,k, [,tau]) 

X 

X  returns  the  absolutely  constant  step  function  with  steps 
X  beta(i) ,  i=l,...,n  and  breaks  0  =  tau(l)  <  ...  <  tau(n+l)  =  1  which 
X  represents  the  (left)-th  coordinate  functional  of  the  B-apline  basis  for  the 
X  knot  sequence  t  :=  [zeros(l.k)  intt  ones(l,k)]  ,  hence  provides  the  norm  of 
X  that  functional  wrto  the  i-norm  on  [0,1]. 

X 

X  On  input,  tau  is  assumed  to  contain  a  reasonable  first  guess. 

X  If  missing  or  inappropriate,  the  knot  averages  are  used. 

tol=l . e~4 ; 

t=[zeros(l,k)  intt  ones(i.k)]; 
npk=length(t) ;n=npk-k; 

X  here  I  omitted  statements  which  would  initialize  tau  as  the  average  of 
X  k  neighboring  knots  in  case  the  initial  guess  provided  is  Inadequate. 

dt=k*diag(ones(l,n) ,/(t(l+k:n+k)-t(l:n))) ;  X  matrix  needed  in  the  computation 

X  of  change  in  tau  . 

b=zeros(n, 1) ;b(left)=k;b(left-l)=-k;  X  generate  the  right  side  ... 
eps=ones(n, 1) ;  X...  and  a  properly  alternating 

eps(l:2:n)=-eps(l:2:n) ;  X...  sequence. 

ttau=tau (2 : n) ; gap=l ; 

iter=0; 
while  (0==0) ; 

X  generate  the  coefficient  matrix.  (Here,  and  below,  spcol(s ,k,tau)  is  the 
X  matrix  whose  i-th  row  consists  of  the  values  at  tau(i)  of  all  the 
X  B-splines  of  order  k  for  the  knot  sequence  s  ,  and  diff(B)  is  the 
X  matrix  with  entries  B(i+1, J)-B(i, J)  .) 

A  =  (diff  (spcoKtt  1] ,k+l,tau))) ' ; 

beta  =  A\b;  X  compute  the  solution  of  the  equation  A*beta  *  b 

betamin=min(abs(beta)) ;  X  compute  the  relative  nonconstancy. . 

gap  =  (max(abs(beta))-betamin)/betamin;  X...  of  abs(beta)  and  ... 
beta-gap  =  [beta(l) ,gap*l.a+4]  X- • •  print  it  out,  along  with  beta(l) 

if  (iter>0)*((gap<tol) I (iter>10)) ,  return;  end 

X  generate  the  change  in  tau; 
c  =  [spcol(t,k,ttau)»dt  zeros (n-1, 1)] ; 

y  =  [-diff(c')  A*eps]\b; 

dtau  =  -(y(l;n-l) ./diff (beta)) ' ; 

ttau=ttau+dtau ; 

X  prevent  changed  tau  from  violating  the  expected  interlacing  by  pulling 
X  back  on  the  Newton  step  if  necessary: 
count-0 ; 

while  (any(ttau<t(l:n-l)) |any(ttau>t(k+2:n+k)) |any(diff (ttau)<=0)) , 
dtau=dtau/2 ; ttau=ttau-dtau ; 

count=count+l;if  (count>20) .error ('no  convergence') , end 

end 

tau=[0  ttau  1]; 
iter=iter+l; 

and 


Figures  5  and  6  parallel  Figures  2  and  3  and  illustrate  thereby  that  the  function  being 
maximized  in  order  to  obtain  Dk  does  appear  to  be  taking  on  that  maximum  at  a  ‘middle’ 
vertex  of  the  domain.  Extensive  calculations  with  the  above  script  for  k  <  21  have  not 
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(^fc+i 1 1 k+2)  £  (0, 1)2 
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Figure  6.  j|  A *  — 1  ,t ||  for  k  =  5  as  a  function  of  the  two  interior  knots 
( tk+l,tk+2 )  €  (0,  l)2 

produced  any  counterexample  to  the  conjecture  that  || A^|S||  is  maximized  when  s  has  no 
interior  knots. 

It  is  also  evident  that  ||  Afc(Sjj  is  minimized  when  s  has  just  one  interior  knot,  of  maximal 
multiplicity,  i.e.,  of  multiplicity  k  —  1.  The  characterization  of  ||A||  as  the  max-norm  of 
A’s  unique  representer  h  in  the  form  of  an  absolutely  constant  step  function  with  n  steps 


12 


makes  it  easy  to  see  that,  in  that  case,  the  norm  is  independent  of  the  location  of  that 
interior  knot. 

Finally,  the  calculations  of  the  representing  step  function  h  presented  no  numerical 
difficulties  in  all  cases  tried  (up  to  k  =  21),  in  marked  contrast  to  the  calculation  of  the 
Chebyshev  splines. 


Conclusion 

There  is  numerical  evidence  that,  in  calculations  devoted  to  bounding  the  p-norm 
condition  number  of  the  (appropriately  scaled)  B-spline  basis,  the  extreme  case  occurs 
for  a  knot  sequence  without  interior  knots,  while  simple  numerical  examples  show  this 
not  to  be  the  case  for  the  max-norm  condition  number  itself.  This  is  disappointing  since 
it  is  only  in  the  latter  case  that  there  seems  to  be  a  formula  available  for  the  condition 
number  when  there  are  no  interior  knots.  Hence,  even  if  the  worst-case  conjecture  for 
the  bound  calculations  for  the  p-norm  condition  number  were  proved,  it  would,  offhand, 
not  help  in  settling  the  problem  of  interest.  This  is  the  proof  that  ail  of  these  numbers, 
Dk,ooi  Dk,i  =  A'(fc),  and  Dk ,  grow  exactly  like  2*,  as  is  suggested  by  numerical  experiment. 
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